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Abstract 

Concepts like peculiar velocity, gravitational force, and power spectrum and 
their interrelationships are of utmost importance in the theories of structure forma- 
tion. The observational implementation of these concepts is usually based on the 
Newtonian hydrodynamic equations, but used up to scales where general relativistic 
effects come in. Using a perturbation of FRW metric in harmonic gauge, we show 
that the relativistic effects reduce to light cone effects including the expansion of 
the universe. 

Within the Newtonian gravitation, the linear perturbation theory of large scale 
structure formation predicts the peculiar velocity field to be directly proportional 
to gravitational force due to the matter distribution. The corresponding relation 
between peculiar velocity field and density contrast has been given by Peebles. Us- 
ing the general relativistic perturbation we have developed, this familiar relation is 
modified by doing the calculation on the light cone in contrast to the usual pro- 
cedure of taking a space-like slice defined at a definite time. The velocity and 
density-spectrum are compared to the familiar Newtonian expressions. In particu- 
lar, the relativistic /3-value obtained is reduced and leads to an increased bias factor 
or a decreased expected amount of the dark matter in a cluster. 



1 Introduction 

Although cosmology has been one of the first areas of application of general relativity, 
in physical cosmology we encounter rarely general relativistic considerations. Relativistic 
effects are commonly restricted to either areas of very high density, like black holes or 
neutron stars, or global cosmological effects, like global dynamics of FRW models. In 
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theoretical implementations of observational cosmology, such as theoretical interpretation 
of the Two-degree Field(2dF)- and the Sloan Digital Sky Survey (SDSS)-projects, the 
Stromlo-APM red-shift survey PJ and the Las Campanas red-shift survey M one usually 
rely on basic Newtonian dynamics [[|, ||, and 0. The correct theoretical interpreta- 
tion of these surveys needs a through understanding of structure formation processes in 
the universe, mainly objects of scales greater than galaxies and at cosmological distances. 
It is obvious that at such scales one can not use the full formalism of general relativity 
to study the dynamics of structures and one has to rely on some approximation methods. 
In relativistic jargon, the Newtonian limit means gravity in the vicinity of objects having 
'small' mass, or expressing it more exactly, where the ratio of the Schwarzschild radius 
to object radius is a very small number 0. But one never discuss the relativistic effects 
on cosmological parameters in the intermediate scales where the relativistic effects due to 
the mass concentration is negligible. This is exactly the case in the theory of structure 
formation. For small density contrasts it is used to apply the Newtonian dynamics to the 
linear theory of structure formation. But even if one tries to apply the linear theory of 
structure formation to evaluate the observational data, general relativistic effects due to 
the large extension of the objects or their distances come in. This is exactly the point 
which has been neglected in theoretical studies up to now and we are going to consider 
it. 



The problem is to formulate an approximation method within general relativity for the 
cases where we have small density perturbation in a FRW model for intermediate regions 
smaller than the Hubble radius but large enough to be forced to consider light cone 
effects. Specifically, we will calculate the relation between the density contrast and the 
peculiar velocity in clusters within the general relativity and compare it to the well known 
Newtonian relations. The formalism developed here shed some light on the methodology 
of the calculation of different so-called light cone effects which has been published |§ , , 
Til , HH > mi- The relativistic calculation helps us to understand more systematically 



the complex role of different cosmological parameters which come in once. In Section 2 we 
review the basic hydrodynamical equations based on the Newtonian gravity. In Section 
3 a perturbation formalism for FRW universes in harmonic gauge is given. Using this 
approximation formalism, we obtain a relativistic relation between the density contrast 
and the peculiar velocity. In Section 4 we obtain relativistic derivation from density 
contrast. Section 5 is devoted to the power spectrum. The general relativistic (3— value is 
then calculated in Section 6. And finally in Section 7 we summarize the main points of 
this study. 

Throughout the paper we choose the signature (—,+,+,-1-), and assume c = 1, except 
where indicated. 



2 Peculiar Velocity Field in the Newtonian Hydro- 
dynamics 

Basic hydrodynamical equations used to study the peculiar velocity of structures in the 
linear regime are the continuity and the Poisson equation: 

5 + V-v = 0, (1) 
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V • g = AnGpS, 



(2) 



where 5 = is the density contrast in the homogeneous and isotropic universe with 
the background density p . The peculiar velocity is defined by v(r) = u(r) — H r, where 
u is the velocity of cosmic fluid, H is the Hubble constant, and r is the physical distance. 
Within the Newtonian gravitational theory, the following relation between the density 
contrast and the peculiar velocity is then derived by Peebles ||: 

where b is the biasing factor defined by 

^galaxy 

and f(tt) = ^2 ~ fi 6 . Eq.(||) can also be written as follows: 

V.»=-« f . (4) 



Note that in general /(fi) depends also on A. However it has been shown by Lahav jTj| 
that the effect of A is negligible. 



3 Peculiar Velocity in the Relativistic Hydrodynam- 
ics 

Our aim is to study the peculiar velocity of a linear structures in FRW space-time. The 
structure is assumed to be large enough to take part in the cosmological expansion, yet 
smaller than the present Hubble radius, with a density contrast small enough to justify 
the use of the linear theory. 

This is in contrast to the familiar textbook approach which uses Newtonian gravity 
for the linear regime in the scales less than the Hubble radius. As a result of our general 
relativistic calculations we will distinguish between intermediate scales where relativistic 
effects in some reasonable approximations have to be taken into account and small scales 
where the Newtonian approximation is valid within the observational accuracy. Therefore, 
we perturb linearly a FRW metric, or what is the same, we expand a general metric linearly 
around a FRW space-time. 

Over the past decades, perturbation of FRW space-times have been studied using different 



gauges like synchronous, Poisson and restricted-gauge (for review see ||14|| , fl"5|j and [[16]). 
These gauges are, however, not suitable for our purposes. As we know from the comparison 
of the weak limit of general relativity to Newtonian gravity, the harmonic gauge is the 
most suitable one playing the role of the Lorentz gauge in electromagnetism as contrasted 
to the Coulomb gauge. Therefore, we try to formulate the corresponding gauge in the 
perturbed FRW space-time. 

Consider a small perturbation of FRW metric in the form 

T„ v = TfJ + 8T^ (5) 
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where gj® is the background FRW metric, Tf$ the stress tensor of the cosmic fluid, h^ v « 
gffi the metric perturbation, and 5T^ V < Tf$ energy momentum tensor of the structures. 
The structures we are considering are in a matter dominated universe. Therefore, the only 
non-zero component of the energy momentum tensor in the comoving reference frame is 
Tqq = Po- Taking into account the general form of the energy-momentum tensor of a 
homogeneous perfect fluid: 

r = (p + /))nV+pf, (6) 

where w M = is the four-velocity of the perturbed cosmic fluid, we may write its 
perturbed component as follows: 

Too = (po + Sp)=p (l + S) (7) 
T 0i = Pou' (8) 
Tij = 0, (9) 

where we ignored higher order of perturbation. We introduce Harmonic gauge in FRW 
metric as follows: 

\K% = V; M (10) 



For the v — and v = i in the expression (10), harmonic gauge can be written in the 
FRW as follows: 

iU°o,o + h\i - h\ + 3-h% - -h\ = v = (11) 
2 2 a a 

3 V» - \v h% - ~h\i + hUj - h° it0 = u = i (12) 

we obtain the (OO)-component of the perturbed Einstein equation 5G° = 8irG5T° , after 
some lengthly calculations, in the following form: 

K w + S-h \ - h w / - 2-h\ - 2{-fh\ - \h l m -3-h\ - l-h 00 , - k-) 2 h 00 = AirGp5 
a ' 2 a a 2 a 2a 2a 

(13) 

Substituting v = component of harmonic gauge expression ([□]) into fljjD, Einstein 
equation obtain as follows : 

- ^oo,/ + 3- Vo - n~ h \o + ?(-) 2/i oo + -(3/ioo - h\) = AixGdp (14) 

Note that for a = constant, i. e. the Minkowski space time, the familiar Poison equation 
is recovered. Now, for scales smaller than the Hubble radius, i.e. A < H~ l = it is 



easily seen that all the terms in the left hand side of (14) can be ignored relative to the 
first term. In fact, we have to compare the Hubble time = " with the characteristic 
time of changes in the cluster, r, which is less that A. It is then obvious that the second 
, third and fourth term which are of the order of {tnT)~ l and fifth term is in the order of 
t~^ are ignorable relative to the first term which is of the order of r~ 2 . Therefore, taking 
the definition of the gravitational potential of the perturbing field, (p, as /too = —20, we 
obtain finally for the perturbed field equation: 

□0 = 4vrGp5, (15) 
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where spatial derivatives of 4> are in the terms of the physical coordinates r. We see that 
the dynamics of perturbing potential obeys the same equation as in the case of relativistic 
radiation, taking into account the difference between comoving and physical-length. We 
may write therefore the potential in the retarded form: 



Gp(t 



)S(t- 



-d 3 r'. 



(16) 



Now, the gravitational acceleration g 
given by: 



-V0 due to the above gravitational potential is 



G(r'-r) . 
rP(t 



G(r' - r) 



)S(t 







2 


c 







+ S(t 



[S(t- 

)p(t 



)p(t- 

}d 3 r'. 



)d 3 r' 



(17) 



Note that the background density in the FRW universe is p oc a" 3 , and in the linear 
regime for the structure formation, we have S oc a. It is then easily seen that the ratio 
of second term in ( ]T7| ) with respect to the first one is of the order of -j^. Therefore, the 
second term may be ignored and we obtain finally 



g(r) 



G(P - r, 



P S(t 



)d 3 r' 



(18) 



One may therefore interpret the gravitational effect of such linear perturbations at inter- 
mediate scales as the Newtonian one with two modifications: 

i. All coordinate lengths should be understood as physical and therefore have to be mul- 
tiplied by the cosmological scale factor a(t). 

ii. The time coordinate should be replaced by the retarded time which means taking into 
account the finite velocity of light or doing the calculations on the light cone. 

This justifies the recent modifications to the Newtonian calculations of different cosmo- 
logical parameters related to the structure formation which are coined with the term 'light 
cone effects', see @, §, 0. 

The expression (|17]) may also be written in the differential form. By taking the divergence 
of (|I7p and ignoring terms of the order of 



s we obtain: 



V • g = 4vrGp5, 



(19) 



which is of the same form as the Newtonian Poisson equation except for the modifications 
discussed above. 

The corresponding continuity equation can be obtained through the Bianchi identities. We 
therefore consider the conservation of energy-momentum tensor in the harmonic gauge. 
Now, using ([F|-|9|), the zero component of the Bianchi identities as the conservation of the 
energy-momentum, 



0. 



(20) 
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leads to 

5 + V-u-^Vo + /4 = 0. (21) 
The third and fourth terms are small relative to the first and second term. This can best 



be seen by rewriting expression ([H]) in the Fourier space: 

6 k + *'Uk- |4§o + *h™ = 0. (22) 

In the linear regime, the first term is proportional to ^5k- From expression (^) it is 
easily seen that the second term is of the same order of magnitude as the first one. To 
estimate the third term we note first that is essentially the potential energy obeying 
the Newtonian relation 

V 2 = AnGp5. 

Taking the time derivative of the corresponding equation in Fourier space we easily obtain: 

4o!o = 3(-) 2 A 2 4-3(^) 3 A 2 4 (23) 
a a 

The two terms on the right hand side of this equation are of the order of (j^) 2 H5k- 
Therefore, the third term in Eq.(|2~l"D is of the order of (j^) 2 times the first or the second 
term and is therefore ignorable. 

To estimate the forth term in (|2lj ) we consider the iO component of the Einstein equation 
by taking into account that A < dn- 

-\h m ,% = 8nGp(^h i0 -u i ) (24) 

Taking the Fourier transformation of the above equation and using the expression (f|), we 
obtain: 

— ^ ~ (—) 2 m k , (25) 

a cLh 

which is small relative to the first and second term and can be ignored. Therefore ([2lD 
may be written as: 

5 + V ■ u = 0. (26) 

We have therefore seen that by perturbing the matter dominated FRW universe in har- 
monic gauge the same Newtonian hydrodynamic Eqs .([[]) and (0) are obtained with the 
two modifications discussed above. Combining the ( |I8| ) and ( p6|) and using the fact that 
in the linear regime 5 oc a we obtain as follows: 

v -"M = Hf v -*M- (2?) 

Omitting V from both sides and substituting g from (0), we obtain for the peculiar 
velocity 

H(t)f(tt) f \r-r'\^ s/ _j \r-r'\^(? -r) 3 , 



If (r, t) = - K ' \ f / Gp(t - -)5(r' } t - 1 U s r' . (28) 

v ' ; AnGbpit) J HK c 1 y ' c '\r-ff 

In the following sections we will rewrite the above relativistic expression in terms of 
observational parameters and will compare it to the corresponding Newtonian results. 
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4 Relativistic Derivation of the Peculiar Velocity From 
Density Contrast 

Assume a cluster of galaxies, G\,G2, ...G n , at a large distance from us. Take two typical 
galaxies G\ and G2 (Fig.l). For the purpose of calculating the relation between density 
contrast and peculiar velocity we need the gravitational action of G2 on G%. The situation 
is best visualized in a space-time diagram. The world lines of G\ and G2 cross our past 
light-cone at the observation event O at P and Q respectively. These events correspond 
to our observation of the galaxies. But the action of G2 on G\ is defined through another 
event: crossing of the past light-cone of P and the world-line of G 2 . Call this event R 
(Fig.l). Times corresponding to each of the events P, Q, and R are characterized by the 
corresponding subscripts. Our observation time at the event O is given by t = to- Now, 
we have to calculate the relation between v pec and 5(x) using observed data on P and Q. 
Note first the following relations between time and space coordinates of the events P and 
Q, taking again c = 1: 

t P = t -\rZ\, (29) 
t Q =t -\r%\. (30) 

Let G\ be the reference galaxy. The action of any other galaxy, like G2, on G\ is along 
the past light-cone of events on the world line of G\. It is easily seen that the time G2 
acts on Gi is given by 

tR = t - \rp\ - \tr - rp\ , (31) 

where rp is defined as the space coordinate of the event R. Within the approximation 
we are interested in, one can replace rp by rq ~ rp. Therefore the above relation can be 
written as 

tR = t - \rp\ - \rg - fp\ . (32) 
Now, Eq.(^) can be written in terms of space time coordinates at P and R: 

*(r P , t P ) = f tr&AUW*. tx) ( ^4k- (33) 

4iiGbo{t P ) J W R - r P \ 



Quantities depending on R have to be reformulated as functions of space time coordinates 
of Q. This is done in two steps. First we expand H(tp) in terms of Ho and p(tp) in the 
terms of p(tp): 

H(tp)=H (l + ^ I ) (34) 



p(tR)=p(t P )(l + 2 lrR TT _ r i pl ), (35) 



where we have used the FRW equations 



a 2 k 8nG 

- + — = -itp 36 

a cr 3 
p + 3H(p + P) = (37) 
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Substituting now the Eqs. 



and (^) into (|33|) , we obtain: 

Hofm f (rfl - rp) ^ 

o{r R ,t R )— —^d r R 



+ 
+ 



H \rp\H f(n) 

c 4irb 
3H H f(Q) 
c Anb 



tr - r P 
I 



\r R - r P \ 



S(r R ,t R ) ;J r S 2 d 3 r R . 
\r R - r P \ 



(3* 



Note that the range of integration is taken as large as desired. Now, terms depending on 
the event R have to be replaced by those depending on Q. Denoting the relative velocity 
of G2 with respect to us as V, we obtain 



f~R = rq - V(t Q -t R ), 



(39) 



where 

V = H f$ + v^ c -v%. (40) 

Here vq is the peculiar velocity of us with respect to CMBR and v pec is the peculiar velocity 
of G2. Using the relations (p0|) and ( p2|) and (0), the expression ( |39"D may be written in 
the form: 

r R = r Q + (v - v - H r Q )(\r P \ + \r P - r Q \ - \r Q \). (41) 

For the large scales we are considering, the term proportional to vq — v pec can be neglected 
relative to HqTq. We may also use the following relation which is easily understood: 



1 



1 



\r R - r. 



\ r Q ~ r P 



l + 3H \r Q \(l--^- 



(42) 



Now, let us expand the density contrast S(t R ,f R ) around (tQ,rQ): 



S(t R , r R ) = 5(t Q , r$ - d5itQ dt ^\ t Q - t R ) + VS(t Q , r® ■ (ft - r$ + ... (43) 

Assuming the density contrast to be proportional to the scale factor a(t), it is seen that 
the time derivation of the density contrast is given by: 



dS{t Q ,rQ) = 2 5(t Q ,rQ) 
~ 3 t Q 



Substituting 



dt Q 

into (f43|), we obtain: 

2 5(t Q ,rQ 



(44) 



S(tR,fh) = 5(t Q ,fQ)-- f 



■{tQ-t R ) + VS-{f R -f Q ). 



(45) 



Substituting for t R and £q from Eqs.flSTf) and ( p2|) , we obtain finally 



5{t Q y Q ) - 2 ^-?\ \rp\ + \r P ~ r Q \ - \r Q \) - V5(t Q ,r Q ) 



3# 

x (|r P | + |r P - r Q \ - \r Q \) - ?Q ' V ^ Tq) (\r P \ + \r P - r Q \ - |r Q |).(46) 



We have still to write the measure of the integral in (|33|), d 3 r R , in terms of the observed 
volume d?TQ\ 

dr R 



d 6 r R 



d 3 r Q . 



(47) 



From Eq.(P), ignoring higher order terms, the Jacobian is obtained to be: 
dr R 



dr. 



Q 



= 1 - V • v(\r P \ + \r P - r Q \ - \r Q \) + v ■ k 
- 3H (\r P \ + \r P - r Q \ - \r Q \) + H f Q ■ k, 



where 



Q 



' Q X P 



(46 



(49) 



\rq-x P \' 

It can easily be seen that, within the approximation we are considering, it is allowed to 
replace in the above relation the divergence of the velocity by the density contrast from 
expression (f|) to obtain the final result: 



d 3 r R 



T$(rQ,t Q )(\r P \ + \rp-T Q \ - \r Q \) - If ■ k 
\r Q \)- H r Q ■ k)d 3 r Q . 



bH~ 

3H (\r P \ + \r P - r Q 



(50) 



Replacing now all terms corresponding to the point R with those defined at the point Q, 
we obtain finally for the peculiar velocity 



v{r P 



H f(Q) f 6(t Q ,r Q ) 



\ r P ~ r Q\ 



[f Q -fp)d 3 r Q + G(r P ), 



(51) 



where 



G(r P ) = F^rp) + F 2 {r P ) + F 3 {r P ) + F 4 {r P ) + F 5 {r P ) 
+ F 6 (r P ) + F 7 {r P ) + F 8 (r P ) + F 9 (r P ), 



(52) 



and 



F x {r P ) 



H 2 f(ty 
Aitb 
r P ■ f Q 

Vq - r p\ 



S{r Q ,t 



- r P 



Q) 



\ r Q - r P\ 



{\rp\ + \r Q -r P \ - \r Q \)(- 



\ r Q\ 



- r P 



F 2 (r P ) 

F 3 (rp) 
F 4 (r P ) 
F 5 (r P ) 



r P 



Airb 



Q) 



^- rP [d 3 
\?Q - rp\ 



Arrb 



*(rq,t Q ) , ( ? ?V r Q , 
\ r Q ~ r P\ 



llHif({l) r 5(t Q ,r Q ) 



3nb 



\rp ~ r Q \ 



(rg - Xp)(\r P \ + \tq - r P \ - \r Q \)d 3 r Q , 



(53) 
(54) 
(55) 
(56) 



Him) ^r^.V ff (^ (r ._ r . )(|rp| + Irg _ rp| _ IrgI)<f (57) 



Anb 



\rp - r Q \ 
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Mr P ) = -^^/^^(M + ir Q -rp|-|r Q |)rgd 3 r Q) (58) 



\ r P ~ r Q\ 



H f(Q) f S(t Q ,r Q ) 
\rp ~ r Q 
H f(Q) r v(r Q ) ■ k 



Fr(r P ) = ^^^4?^ / T^^^^^T^ ^ I I ~^ ^ ' ^ | ^ | ) ^ C ^5 ) , (59) 

47T0 J rp-rn 



Fs{rp) = + 4ri I ■- - ^ Q-^ntQ^Q)d ^ (60) 



\ r P ~ tq\ 



H m) r v(r Q )-V6(t Q ,r Q ) 
Airb J \r P — tq| 3 

To simplify the results and to have an estimation of the corrections to the Newtonian 



- H f(n) f y(r Q ) ■ V^(t Q ,r Q ) _ ^ 3 
F ^ rp > = ArtT J \r P -r Q \ 3 (rQ-rp){\rp\ + \r Q -rp\-\r Q \)d r Q (61) 



expressions, we take a closer look at different terms in fl5"2|). Let us first distinguish 

between F\ — F±, and Fq on one side and F-j — Fg on the other side. The latter terms 

are typically of the order of formers times the fraction of the peculiar velocity divided 

by the Hubble velocity, being of the order of 1CT 2 . Therefore, the latter terms may be 

neglected. The term F§ contains the factor tq ■ V5 which is vanishing in general. To 

see this, consider a cluster around Tq. Because of the symmetric distribution of matter 

in the cluster the gradient term accept positive and negative values while rg is almost 

constant. Therefore, one may assume that the integral in F 5 is vanishing. The same is 

true about the contribution of the term rQ ' rp - 2 in Fx- This is due to the large extent of 

|fp-rb| 

our integration domain represented by rq. Now, for Stromlo-APM- and Las Campanas- 
red-shift survey we may assume the distance to the cluster to be of the same order as the 
extent of the cluster, L, so that we have \ fp — Tq\ ~ \ fp\ ~ |fg| ~ L ~ lOOMpc. A closer 
look at the remaining terms shows that Fx, 3,4,6 are of the order of or bigger than F2. We 
may therefore write 

G ^ ^W, (62) 

where 

7 H f(n) f S(t Q ,r Q ) _ 3 

47T0 J \rp — Vq\ 

Substituting now G in expression (|5l"D we obtain finally: 

v^W(l + ^-), (64) 

where L is the size of structure and c is explicitly inserted again. In the limit c — > 00 or 
L — > we obtain the Newtonian value for the peculiar velocity, vpj, which is just the first 



term on the right hand side of (51). Therefore, the relative relativistic correction to the 

(65) 



peculiar velocity up to the first order of - is given by 

Vrel -VN 4 LH 



vn 3 c 

where we have taken the absolute values for the velocities and added the subscripts re ; 
and 7v to emphasize the relativistic corrections. Taking the divergence of the ([51]), we 
obtain a relation between the density contrast and the divergence of the peculiar velocity: 

^V-ffM^M-^V-GM. (66) 



10 



Here again, in the limit c — > oo the second term on the right hand side vanishes and we 
get the familiar Newtonian expression for the density contrast, which will be denoted by 
<5jv( r p)- Now, to distinguish the density contrast appearing on the right hand side of 



from the Newtonian value, we call it 5 re i. The expression (|66|) is now written as: 



Srei(r P ) = 5 N (r P ) + & V ■ G(r P ). (67) 
H f{U)c 

In the following sections we will calculate the density and velocity power spectra for both 
Newtonian and relativistic cases. 



5 Power Spectrum: Newtonian Versus Relativistic 

Case 

Take the correlation function of density contrast as 

«r) = Js(r + r')S(r')^, (68) 

defined in a cosmological volume V. For simplicity we will omit hereafter the subscript p 
in rp. The power spectrum is given by the Fourier transform of the correlation function: 

f (r) = J P{k)e tkr d 3 k. (69) 

Using relation fl6"T|) between the Newtonian and relativistic density contrast, substituting 
it in Eq.flBlf) and ignoring higher order terms in 5, we obtain the following relation between 
the Newtonian- and relativistic-correlation functions: 

6r(r) =M')-^4ny / 8(r ' + r)\ >-G(r')d 3 r ' -^j^j J 5{r')V -G{r' + r)dh' . (70) 

The divergence term on the right hand side is calculated in the Appendix A and is given 
by 

V-G{r)= ~ AH ll m 8{r,t) \r\ (71) 

Using this relation and taking the Fourier transform of expression (|70|) , the following 
relation between the Newtonian- and relativistic-power spectrum is easily obtained: 

P N (k) = P(k) rel (l + ^^)- (72) 
In the limit c —>■ oo, the Newtonian value is regained. 

6 (3— Value: Newtonian Versus Relativistic 

Let us now calculate the relativistic velocity power spectrum. The peculiar velocity is 
now expressed in terms of its Fourier components: 

v(r) = ^ / v k e^d% (73) 
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where V is the space volume. Its power spectrum is defined as Pv(k) =< \vk x \ 2 + \vk y \ 2 + 
\vk z \ 2 >■ If v(r) is an isotropic Gaussian field, then the different Fourier components are 
uncorrelated and the power spectrum provides a complete statistical description of the 
field. The Velocity spectrum is given by 

V\k) = ^k 3 P v (k). (74) 



As we have shown in the appendix B, the Fourier transform of (|65| ) leads then to the 
following relation between the Newtonian- and relativistic- velocity spectrum: 

= V&(*)(1 + ^). (75) 

The Newtonian spectrum is independently obtained by taking the Fourier transform of 
(I): 

V^(k) = ^- 2 P 2 HlkP(k). (76) 

where (3 = ^p-- These relations are obtained in terms of real space expressions in contrast 
to power spectrum data which are obtained in the red-shift space |17|]- It is used to change 
the factor f3 2 in front of the term on the right hand side of ( |7E|) to F(/3) to account for 
the transformation between red-shift- and real-space. Therefore we may write the final 
equation in terms of red-shift data in the form 

Vm = ^- 2 FU(3)HokP(k), (77) 

where 

F »<"> - i + wf s+ w (78) 

Now, using expression ( [75] ) the relativistic velocity spectrum will be given in the form 

V 2 e i{k) = ^- 2 F^ el ((3)H 2 kP(k), (79) 



where 



F 2 M = FU(3)(l+ 1 -^)- (80) 



7 Discussion 



Any observation in cosmology is carried along the light cone and not on a space-like slice 
defined by a definite observer time. However the difference between the two procedures 
are usually ignorable, which leads to the usual ignorance of the light cone effects and 
the finite signal velocity. At cosmological distances this difference can leads however to 
observable effects. 

Here we have calculated the peculiar velocity-density contrast relation in a cluster of 
galaxy taking into account the finite signal velocity. The result differs from the familiar 
Newtonian relation through a lengthly term which can however be simplified for a typical 
cluster. The estimation based on this simplification depends on the extend of the cluster; 
the larger the extent the bigger the difference. 
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Appendix A 



Let us calculate the term V • G(r) appearing in ( |67|) and (ff5|) . Using eqs. (|53| - |6TD we may 
write: 

ggg/gg /• rp-r q 3 3ff 2 /(^) /• |r P | 

4:7Tb J \rp\\rp — tq\ 4tto J \rp — tq\ 

_ «f°> M ,( rp ), (8!) 

The symmetric distribution of matter in large scales leads to positive and negative values 
for rp ■ tq in the first integral which makes it vanishing. The second integral on the right 

hand side may be written as 3H " ^ n ~ > \rp\ J S( ~ rQ \a d 3 r' ~ 3H ° < 5 > L. The average 

\rp — rQ 47to 

value of density contrast, < 5 >, taking the domain of integration sufficiently large, does 
vanish too. We then obtain finally: 

V.d(r P ) = \rp\6(r P ). (82) 

3b 

Appendix B 

In order to derive the Newtonian power spectrum we substitute expression (|7lf) into (|70[) : 

6r(r) = 6w(r) + ^ / <Kr' + r)5(r') \r'\ d 3 r f . (83) 



3c 

Taking the Fourier transformation of it we obtain 
J P N {k)e lk - r d 3 k = J P Rel {k)e lk - r d 3 k 



+ j^y 3 ^ / 5 k 5' k e lk -^e lk '- r 'r'd 3 r'd 3 kd 3 k'. (84) 

Putting for simplicity r' = ^ we may write / e ir ■( k + k ') ( j l s r l = (27r) 3 5 3 (fc + k'), which leads 
to 

P N (k) = P Rel (k)(l + ^^). (85) 

Calculation of the velocity power spectrum is similar to above procedure. One start with 
(|65|), by substituting it into (]73D and taking Fourier transformation of it and using Eq.fl74|) 
, it is easily seen that: 

^)=^)(1 + ^). (86) 
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8 Figure Captions 

Fig. 1. Our past light cone crosses the world-lines of galaxies G\ and Gi at P and Q 
respectively. The event R, on the intersection of the world line of G2 and the past light 
cone of P, gives the moment at which G2 acts on G\ being effective at P. 
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